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Abstract. We numerically estimate the semi-major axis chaotic diffusion of the Vesta family asteroids induced 

by close encounters with 11 massive main belt asteroids : (1) Ceres, (2) Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) 
Hygiea, (15) Eunomia, (19) Fortuna, (324) Bamberga, (532) Herculina, (704) Interamnia. We find that most of 
the diffusion is due to Ceres and Vesta. By extrapolating our results, we are aide to constrain the global effect 
of close encounters with all the main belt asteroids. A comparison of this drift estimate with the one expected 
for the Yarkovsky effect shows that for asteroids whose diameter is larger than about 40 km, close encounters 
dominate the Yarkovsky effect. Overall, our findings confirm the standard scenario for the history of the Vesta 
family. 



1. Introduction 

It is now admitted that most of V-type near-Earth aster- 
oids (NEAs) and howardite, eucrite and diogenite (HED) 
meteorites are fragments of a collision between Vesta and 
another object which also produced the Vesta family (see 
Sears 1997). V-type NEAs and HED meteorites are former 
members of the Vesta family that have been carried from 
the main asteroid belt to Earth-crossing orbits by the two 
main resonances in the neighborhood of Vesta : the vq and 
3 : 1 resonances with Jupiter. The major issue in this sce- 
nario is the fact that the average life time of fragments 
in these two resonances is too short to explain the mean 
age of HED meteorites and V-type NEAs (Migliorini et al. 
1997). 

The commonly accepted explanation for this is that 
NEAs spent most of their life time in the main asteroid 
belt before entering in resonance and being carried to their 
current orbits. This supposes that a mechanism can induce 
a diffusion process that brings fragments to one (or both) 
of the two main resonances. The Yarkovsky effect is now 
admitted to be the main mechanism explaining the diffu- 
sion of Vesta family members (e.g. Carruba et al. 2003). 

However other processes also contribute to the diffu- 
sion. It is in particular the case for the chaos induced by 
overlaps of mean motion resonances or close encounters 
between asteroids in the main belt. In the case of the res- 
onances, Morbidelli & Gladman (1998) showed that the 
chaotic zone around the 3 : 1 resonance is too narrow to 
explain the needed diffusion. Recently, Laskar et al. (2011) 
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showed that close encounters among massive asteroids in- 
duce strong chaos which appears to be the main limit- 
ing factor for planetary ephemeris on tens of Myr. In the 
continuation of this work we concentrated our study on 
the diffusion induced by close encounters of Vesta family 
members with massive asteroids. 

This effect has already been studied for different as- 
teroid families (e.g. the Flora and Lixiaohua families 
Nesvorny et al. 2002; Novakovic et al. 2010). The most 
exhaustive study (in terms of number of asteroids taken 
into account) concerned the Gefion and Adeona families 
(Carruba et al. 2003). The authors considered all 682 as- 
teroids of radius larger than 50 km for encounters with 
members of the families. They concluded that the four 
largest asteroids ((1) Ceres, (2) Pallas, (4) Vesta and (10) 
Hygiea) had a much larger influence than all the 678 re- 
maining asteroids so the latters can be considered as negli- 
gible. The Vesta family has also been the object of such an 
analysis (Carruba et al. 2007) but only close encounters 
with (4) Vesta were taken into account. These different 
studies show that the elfect of close encounters with mas- 
sive asteroids depends a lot on the considered family and 
its position in phase space with respect to those of massive 
asteroids. 

The purpose of this work is to explore further the 
effect of close encounters for the Vesta family by con- 
sidering the effect of 11 large asteroids : (1) Ceres, (2) 
Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) Hygiea, (15) 
Eunomia, (19) Fortuna, (324) Bamberga, (532) Herculina, 
(704) Interamnia. Moreover, we use the results from these 
interactions to extrapolate and obtain an estimate of the 
effect which would be raised by the whole main asteroid 
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belt. Finally we compare the semi-major axis diffusion in- 
duced by the Yarkovsky effect and the diffusion due to 
close encounters (in the case of the Vesta family). 

2. Numerical simulations 

We ran two sets of numerical simulations. Both of them 
comprise the 8 planets of the Solar System, Pluto, the 
Moon and the 11 asteroids enumerated above. In both sim- 
ulations, we generate a synthetic Vesta family produced by 
an initial collision of another object with Vesta. We do not 
model a realistic collision process but instead we consider 
a set of test particles initially at the same position as Vesta 
but with different relative velocities. This is not a problem 
because we are mainly interested in this work in exploring 
the phase space in the region of the Vesta family. The rel- 
ative velocities of the fragments of collision are sampled 
both in norm and direction. For the norm, we take 8 dif- 
ferent values (every 50 ms~^) between the escape velocity 
at Vesta's surface (around 350 ms~^) and twice this value 
(700 ms~^). The inclination is sampled between —50° and 
50° every 10° (11 different values). Finally, the direction 
in the invariant plane is taken every 10° (36 values). This 
way we create 3168 test particles that represent the Vesta 
family. We plotted the positions of these particles in the 
(a, e) and (a, /) planes (where a is the semi-major axis, e 
the eccentricity and / the inclination) superimposed with 
the catalog published by Nesvorny (2010) of current Vesta 
family members (Fig.l). Note that despite we did not com- 
pute the proper elements whereas the catalog uses them 
the two sets of points superimpose quite well. Moreover we 
constructed an initial set of particles just after the origi- 
nal collision whereas the catalog exposes the remnants of 
the family after approximately 1.2 Gyr of evolution (see 
Carruba et al. 2007). 

In the first simulation (the reference simulation, S), 
the 11 asteroids are considered as test particles as well 
and thus do not perturb the fragments of collision. In the 
second simulation (SE) the 11 asteroids are considered as 
the planets and thereby interact with members of the syn- 
thetic Vesta family. Both solutions are integrated over 30 
Myr using the symplectic integrator described in (Laskar 
et al. 2011) and references therein. The effects of all plan- 
ets and Pluto are taken into account, as well as general 
relativity. 

In order to be able to analyze the evolution of the test 
particles and in particular the impact of close encounters 
with the 11 asteroids we set up different logs. We record 
for each particle its instantaneous orbital elements every 
1 kyr. We do not compute proper elements but instead we 
use the minimum and maximum values (see Laskar 1994) 
taken by a, e and / on time spans of 10 kyr. Actually, we 
use the average of the minimum and maximum values as 
proper elements. We also set up a log of close encounters 
for each particle. Each time that a particle passes within 
0.01 AU from an asteroid, the asteroid number, the mini- 
mum distance of approach and the time of the encounter 
are recorded in the logs. 
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Fig. 1. Initial distribution of the synthetic Vesta family 
in the (a, e) plane (top) and the {a, I) plane (bottom). 
Real members proper elements (from Nesvorny 2010) are 
superimposed for comparison. 

3. Global overview of the diffusion 

The goal of this study is to understand whether the cur- 
rent flux of V-type asteroids coming from the main belt 
and carried to near-Earth orbits through strong planetary 
resonances can be explained by the chaos induced by close 
encounters in the main asteroid belt (or at least that these 
encounters contribute significantly). The simplest way to 
check this assumption is to look in both S and SE simula- 
tions to the number of NEAs as a function of time and to 
see if the consideration of asteroidal interactions implies 
a greater number of NEAs. We use the usual criterion to 
decide whether an asteroid is a NEA or not : if the perias- 
tron of its orbit become smaller than 1.3 AU the asteroid 
is considered as a NEA (see for example Morbidelli et al. 
2002). In Fig. 2 we plot the evolution of the number of 
NEAs in both simulations as a function of time. We can 
see on this graph that the close encounters have not a 
major effect on the population of NEAs on the duration 
of the simulation. However, we did not ran our simula- 
tions on the total age of the Vesta family (around 1.2 Gyr) 
and we can think that just after the collision that origi- 
nated the Vesta family, the population of V-type NEAs is 
dominated by fragments of the collision that were directly 
injected in the two main resonances. The diffusion process 
is supposed to provide these resonances on a much larger 
time scale. 

Another way to look for the effect of close encounters 
is to compare the proper elements of the initial and final 
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Fig. 2. Number of near-Earth asteroids (NEA) among the 
synthetic Vesta family as a function of time in S and SE. 
The criterion for considering a fragment as a NEA is based 
on the minimum value taken by the periastron. When this 
value is less than 1.3 AU the fragment is considered as a 
NEA. 

distribution of fragments in both simulations. A usual way 
to do that is to compute the standard deviation of different 
proper elements (and in particular the semi-major axis) on 
the whole population of fragments at a time t with respect 
to the initial values of these elements (see for example 
Carruba et al. 2003) : 



o"i(i) = 



N-l 



(1) 



where x can be the semi-major axis a, the eccentricity e 
or the inclination / of Vesta family fragments. As we men- 
tioned it we use the averages of the minima and maxima as 
proper elements and in particular for the initial values we 
compute them on the first 10 kyr of the simulation. We do 
not use the initial conditions given in Fig.l because they 
are only instantaneous values. 

We measure the diffusion in semi-major axis only due 
to close encounters, thereby we do not take into account 
resonances and in particular strong ones. For this purpose 
we plotted the value of the standard deviation between the 
beginning and the end (after 30 Myr of evolution) for both 
simulations {S, SE) as a function of the initial (proper) 
semi-major axis. We divided the interval of semi major 
axis in bands of 0.01 AU and for each band we computed 
the standard deviation (of the semi-major axis, the eccen- 
tricity and the inclination) of the set of asteroids that are 
initially in this band. The results of this calculation are 
given in Fig. 3. We can see that the strong resonance 3 : 1 
around 2.5 AU induces an important diffusion in both S 
and SE and that for semi-major axis lower than 2.25 AU 
both simulations are strongly affected by resonances (see 
Nesvorny et al. 2008, for a list of resonances acting on the 
Vesta family). Note that in strong resonances most frag- 
ments are highly unstable and for a great part of them the 
calculation stops before the end of the simulation (colli- 
sions with a planet or the Sun or escapes from the Solar 
System) . Thereby calculations of the diffusion between the 
beginning and the end of the simulation in such zones is 
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Fig. 3. Semi-major axis dependency of the diffusion in 
semi-major axis (top), eccentricity (middle) and inclina- 
tion (bottom) in S and SE. For each band of 0.01 AU, 
we plot the standard deviation of the average of the mini- 
mum and the maximum values of a (respectively e and /) 
during the last 10 kyr of the simulations with respect to 
the initial values (see Eq.l). 



not representative of their instability. For instance the dif- 
fusion is zero between 2.5 and 2.51 AU only because none 
of the fragments that were initially in this zone finished 
the simulation. 

Between 2.26 and 2.48 AU, S shows a very limited 
semi major-axis diffusion even if we can see small peaks 
(e.g. around 2.42 AU corresponding to the 1 : 2 resonance 
with Mars). SE undergoes a more important diffusion in 
this band. We run our calculations on this band because it 
fulfills at the same time two important conditions : the dif- 
fusion due to resonances is limited and the band contains 
(at the beginning of the simulation) a sufficient amount 
of fragments (2617 of the 3168 fragments) which is impor- 
tant for the statistics. Note that for the eccentricity and 
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Fig. 4. Evolution in time of the semi-major axis diffusion 
in S and SE. We plot the standard deviation of the av- 
erage of the minimum and the maximum values of a for 
each fragment over 10 kyr steps with respect to the ini- 
tial values (see Eq. 1) for asteroids whose semi-major axes 
are initially between 2.26 and 2.48 AU (in order to avoid 
strong resonances). 



the inclination, the diffusion is hidden in the remaining 
oscillations due to the secular evolution of these elements 
induced by the planets and that we cannot see any dif- 
ference between both simulations in Fig. 3. However, this 
is not really a problem since in this work we are mainly 
interested in the semi-major axis diffusion. 

The evolution (in time) of the standard deviation of 
the semi-major axis computed on the 2.26 - 2.48 AU band 
for both simulations is plotted in Fig. 4. The standard de- 
viation is approximately constant for S which reinforces 
our choice of this band. On the contrary, SE clearly under- 
goes diffusion. Analyzing Fig. 4 gives us an approximation 
of the drift in semi-major axis due to asteroidal interac- 
tions during the simulation. After 30 Myr of evolution, we 
have a standard deviation around 9 x 10"** AU for SE 
and 2 x 10"'* AU for S. In order to evaluate the diffusion 
due to asteroidal interactions, we have to remove the dif- 
fusion obtained in the reference simulation (which is prob- 
ably due to resonances and output sampling of the secular 
evolution of the semi-major axis) from the total standard 
deviation. Note that the total variance (the squared of the 
standard deviation) is the sum of the variances of the dif- 
ferent contributions. Thus we obtain a standard deviation 
due to asteroidal interactions after 30 Myr of evolution of 
about 8.8 X 10~* AU. We also observe several jumps in 
the semi-major axis diffusion, which are due to very close 
encounters of a single fragment with one of the 1 1 consid- 
ered asteroids. Actually it has been shown (Carruba et al. 
2007) that the distribution of jump sizes is not a Gaussian 
but has thicker wings. Thus the curve is not regular on 
the time scale of this simulation because of the presence 
of rare very important events. We can find the same kind 
of effects on Fig. 3, the curve is not regular in the 2.26 - 
2.48 AU band and shows peaks due to a few very close 
encounters. 



=1 
< 




Fig. 5. Semi-major axis evolution of a fragment which un- 
derwent various close encounters. We plot the maximum, 
minimum and instantaneous values of the semi-major axis 
of the fragment. All times of encounters closer than 10~^ 
AU are highlighted with dotted vertical fines. 
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Fig. 6. Semi-major axis evolution of a fragment whicfi 
underwent a very close encounter (minimum distance of 
3.4 X 10"^ AU) with Vesta (at t w 8.39 Myr). We plot 
the maximum, minimum and instantaneous values of the 
semi-major axis of the fragment. All times of encounters 
closer than 10"'^ AU are highlighted with dotted vertical 
lines. 



4. DifFusion induced by asteroids close encounters 

The global analysis that we made does not allow us to 
separate the contributions of the 11 asteroids in the diffu- 
sion process and to extrapolate to the whole main asteroid 
belt. In order to do that we have to examine logs of close 
encounters and estimate the diffusion resulting from each 
of these events and then make statistics. All this analysis 
is still restricted to the 2.26 - 2.48 AU band in order to 
avoid the strong resonances as previously. 

For each encounter recorded in the logs we compare 
the value of the semi-major axis before and after the en- 
counter. Fig. 5 gives an example of the evolution of the 
semi-major axis of a fragment that underwent several close 
encounters which resulted in jumps of different sizes. Fig. 6 
gives a more unusual example of the evolution of a frag- 
ment which underwent a very close encounter with Vesta 
resulting in a very important jump in semi-major axis. 
Then for each of the 11 asteroids we calculate the stan- 
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dard deviation of the distribution of jump sizes induced by 
close encounters with this asteroid. This gives a measure 
of the average diffusion resulting from a single encounter 
with the asteroid. As we already noticed, it is more con- 
venient to manipulate variances than standard deviations 
since the latter are not additive. If we assume that close 
encounters are independent events, every fragment under- 
goes a random walk and after TV encounters, the total 
variance is multiplied by a factor N : 

allN] = Nal[l]. (2) 

This gives a measure of the average diffusion resulting 
from A'' encounters with the considered asteroid. We just 
have then to replace N by the mean number of close en- 
counters per fragment with the selected asteroid during 
the whole simulation and we obtain the diffusion (vari- 
ance) due to this asteroid on the duration of the simula- 
tion. The diffusion due to close encounters with all the 11 
asteroids on the duration of the simulation is the sum of 
the variances of the 11 asteroids. With the same reasoning 
as for Eq.2 it is possible to extrapolate the value of the 
variance on longer time scales : 

a^it) = -^al{T,im). (3) 

sim 

Note that it as been shown (e.g. Nesvorny et al. 2002) 
that close encounters do not exactly result in a random 
walk (there are some correlations between successive en- 
counters) and that the variance is not exactly linear with 
time (or number of encounters). It can, indeed, be written 
as a power law of the form : 

aa(t) = C7i«. (4) 

For a random walk : B = 0.5. Numerical estimations of 
this coefficient B for different families (see Nesvorny et al. 
2002; Carruba et al. 2003; Novakovic et al. 2010) show 
that it depends on the family considered and can even 
evolve with time. However the values found in the litera- 
ture range between 0.5 and 0.65 and we did not find any 
estimations in the case of the Vesta family so we decided 
to stick to the random walk hypothesis. 

Note that in this reasoning we do not take into account 
the influence of the noise in the calculation of jumps sizes. 
Indeed, we observe on Fig. 5 and Fig. 6 that the minimum 
and maximum values of the semi-major axis are not con- 
stant outside the close encounters (jumps) but undergo os- 
cillations which are due to remaining secular evolution. It 
is possible to eliminate an important part of this noise by 
taking the minimum of minima and the maximum of max- 
ima over greater time scales before and after the encoun- 
ters. In order to choose the best time interval to compute 
these extrema and to evaluate the remaining noise after 
this treatment we use the reference simulation S and we 
calculate the standard deviation of the difference of semi- 
major £ixis before and after random times. Since there 
are no close encounters thus no jumps in this simulation, 
this standard deviation measures only the remaining noise. 



With this method we compute the standard deviation of 
the noise for different time intervals for the calculation of 
extrema and we find that the minimal value (4.4 x 10~^ 
AU) is obtained for an interval length of 200 kyr. This is 
the value we employ for all the following calculations. 

When we compute the variance of jump sizes in SE 
(for real close encounters) we are affected by the same 
noise. The value we obtain for the variance is the sum of 
the variance of the real diffusion and the variance of the 
noise. Thus the real diffusion resulting from close encoun- 
ters is given by : 

^aW^'^aW-^Lae- (5) 

Another problem we experienced is that the logs of 
close encounters record all encounters within 0.01 AU but 

it appears that at this distance the effect is too small and 
is completely hidden in the noise. Thus statistics are con- 
taminated by false events (with no real jumps). In or- 
der to avoid this contamination we take into account only 
encounters within 0.001 AU which are more significant 
events. 

Finally, when a fragment undergoes several close en- 
counters spaced by less than 200 kyr, we are not able to 
separate the effect of each encounter. We thus ignore all 
such multiple encounters and we keep only single encoun- 
ters in the calculation of the standard deviation of the size 
of the jumps. However, when we compute the total diffu- 
sion on the duration of the simulation we use the total 
number of encounters that occur during the simulation as 
multiple encounters should not be ignored. 

Table 1 sums up the results we obtain for each asteroid 
and for the 11 asteroids together. We give the numbers of 
encounter during the simulation (including multiple en- 
counters), the standard deviation obtained from statistics 
on single encounters ((Tail]): the corrected values from the 
noise (ffall]) and the standard deviations on the duration 
of the simulation (aa(Tsj„)). Finally, the percentages of 
contribution of the different asteroids are given in term of 
variance. We obtain a total diffusion due to close encoun- 
ters of 8.63 X 10"* AU during the 30 Myr of the simulation 
(see table 1). This number is to be compared with the drift 
rate obtained with the global approach (8.8 x 10~^ AU). 
The results given by both methods are in good agreement. 

By using Eq.3 we can extrapolate our results to longer 
times : 

aa{t) = 1.57 X IQ-^y/t (Myr) AU. (6) 

Of course it is possible to do the same for the contribution 
of each asteroid. This allows us to compare our results with 
the previous study of Carruba et al. (2007). In this article, 
the authors find a drift rate of 2.0t2:o ^ 10"^ AU/(100 
Myr) by considering only close encounters with Vesta. If 
we extrapolate to 100 Myr and consider only Vesta we find 
a drift of 1.3 x 10~^ AU which is coherent with Carruba 
et al. (2007) findings. 

If we look at the different contributions in table 1 we 
see that close encounters with Ceres and Vesta together 
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Table 1 . Comparison of the contributions of the encoun- 
ters with the 11 considered asteroids to the semi-major 
axis diffusion. All standard deviations are given in AU. 
The contributions are given as percentages of the total 
variance (which is additive unlike the standard deviation) . 
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represent about 99% of the total diffusion (in variance). 
Contributions from the 9 other asteroids seem negligible. 
This is coherent with findings of previous studies on other 
asteroid families (see in particular Carruba et al. 2003). 
On the contrary, since Ceres represents 36% of the diffu- 
sion this asteroid should not be neglected (like in Carruba 
et al. (2007)). We believe that two main reason can ex- 
plain these proportions of contribution : the mass of the 
asteroids and their proximity with the Vesta family mem- 
bers in phase space. The mass must influence the impor- 
tance of the effect of a single close encounter whereas the 
proximity in phase space influences the number of such 
events. In order to check for these correlations we plotted 
the variance due to close encounters as a function of the 
mass of the considered asteroid (Fig. 7) and the number of 
events as a function of the distance in phase space of the 
considered asteroid with respect to Vesta (Fig. 8). We use 
the definition of the distance in phase space introduced by 
Zappala et al. (1994) : 



d: 



+ 2{Sey + 2{5siniY 



(7) 



This distance has the dimensions of a velocity since it is an 
estimation of the ejection velocity that would be needed 
to carry two objects initially at the same position to their 
current positions. In Eq.7, we used the mean values of the 
semi-major axes, the eccentricities and the inclinations of 
the 11 asteroids during the whole simulation (30 Myr) in 
order to compute the distances. We deduce from these 
plots that we indeed have an effect of the mass and the 
proximity in phase space but we cannot obtain a simple 
and precise law in order to evaluate the diffusion that 
would result from close encounters with other asteroids 
than the ones considered in our simulation. 

However, we can constrain the global effect of the re- 
maining objects of the main asteroid belt. In our simula- 




Fig. 7. Contribution to the semi major axis diffusion due 
to each considered asteroid as a function of its mass. 




10 15 

d (AU/yr) 

Fig. 8. Number of close encounters (< 10"'^ AU) between 
Vesta family fragments and each of the 11 considered as- 
teroids as a function of the distance in phase space be- 
tween these asteroids and Vesta. We use the mean values 
of the semi-major axes, the eccentricities and the inclina- 
tions of the 11 asteroids during the 30 Myr of the simula- 
tion in order to compute the distances (see Eq.7). 



tion we took into account the most massive objects of the 
main asteroid belt. Actually, if we add the masses of the 
11 asteroids of our simulation we obtain a total mass of 
8.36 X 10"1°Mq. The estimated total mass contained in 
the main asteroid belt is about 15 x 10~^°Mq (obtained 
from INPOPlOa fits, Fienga et al. 2011) while Krasinsky 
et al. (2002) gave 18 x 1O"^°M0. The most massive as- 
teroid in the main belt is Ceres with 4.76 x IO-^^Mq, 
thus the remaining mass (not considered in our simula- 
tion) represents about 1.5 to 2 times the mass of Ceres. 
The total mass contained in the main asteroid belt is not 
well constrained but we can assume the upper limit of this 
remaining mass to about twice the mass of Ceres. 

In Fig. 7, the slope of the regression curve is about 1.2. 
We plotted a straight line of slope 1 for comparison. A 
slope of 1 means that two objects of masses 0.5 have the 
same effect than one of mass 1. If the slope is higher than 
1 it is more efficient to have a single object whereas if 
the slope is lower than 1 it is more efficient to have two 
objects. Even if this slope is not well constrained in Fig. 7 
it is reasonable to consider that it is greater than 1 so 
it is more efficient to have one big object than to split 
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it into smaller objects. Thereby, we can assume that the 
total contribution of the remaining objects of the main 
belt will be less than twice the contribution of Ceres. 

We still have the problem of the influence of the dis- 
tance in phase space but it seems realistic to consider that 
the case of Vesta is singular since it is the parent body of 
the Vesta family and that the remaining objects must have 
a probability of close encounter closer to Ceres 's one than 
Vesta's one (see Fig. 8). Thus, the total variance due to 
encounters with all the asteroids of the main belt should 
not exceed (and is probably a lot smaller than) 1.7 times 
the value of the variance obtained in our simulation. In 
terms of standard deviations it corresponds to a factor of 
about 1.3. 

5. Comparison with the Yarkovsky- YORP effect 

The Yarkovsky effect is a non-gravitational force which 
results from the interaction between asteroids and the 
solar radiation. There are actually two versions of the 
Yarkovsky eS^ect : a diurnal and a seasonal effect. The 
diurnal effect can result in an increase or a decrease of 
the semi-major axis depending on the obliquity (e) of the 
asteroid (oc cose), whereas the seasonal effect systemati- 
cally decreases the semi-major axis (oc sin^ e). The diurnal 
effect is maximal when e = 0° or 180° and is zero when 
e = 90°. On the contrary, the seasonal effect is maximal 
when e = 90° and zero when e = 0°, 180°. Both effects are 
size dependent, for instance the diurnal effect is given by 
(see in particular Nesvorny et al. 2008) : 



events is given by (see Farinella et al. 1998; Morbidelli & 
Vokrouhlicky 2003) : 



da 



= 2.5 X 10" 



1 km 
D 



cose AU/Myr. 



(8) 



In addition to the Yarkovsky effect, we have to consider 
the YORP (Yarkovsky - O'Keefe - Radzievskii - Paddack) 
effect which also comes from the solar radiation but affects 
the rotational velocity and the obliquity of asteroids. This 
effect depends a lot on the shape and the surface compo- 
sition of the asteroid but it is possible to do some statis- 
tics to understand what is the most probable scenario 
(see Vokrouhlicky & Capek 2002; Capek & Vokrouhlicky 
2004). It has been shown (Capek & Vokrouhlicky 2004) 
that for most of basaltic asteroids (more than 95%), the 
obliquity is asymptotically driven to 0° or 180° by the 
YORP effect. The time scale for the YORP effect is about 
10-50 Myr (see Vokrouhlicky & Capek 2002: Morbidelli & 
Vokrouhlicky 2003). This means that on a 1 Gyr time scale 
(which is the order of magnitude for the age of the Vesta 
family (Carruba et al. 2007)), we can consider that the 
obliquity is either 0° or 180° during the whole evolution 
and that the seasonal effect is negligible and the diurnal 
effect is maximal (in one or the other direction). 

Nevertheless, the collisions must also be taken into ac- 
count. We have to distinguish whether a collision is de- 
structive or not. The collisional life time of a 1-10 km as- 
teroid in the main belt is about 1 Gyr. But non-destructive 
collisions are more frequent and can reorient the spin 
axis of the asteroid. The characteristic time scale of such 



For kilometer-sized asteroids, the time between two reori- 
enting collisions is about 300 Myr, it means that there are 
a few (3 or 4) reorientations on the time scale of the age 
of the Vesta family. Thus it is likely that these reorienta- 
tions will reduce the diffusion due to the Yarkovsky effect 
but not by a large factor. Moreover, we can still consider 
that the YORP effect acts almost instantaneously between 
two reorientations and that the seasonal Yarkovsky effect 
is negligible and the diurnal effect is always maximal. 

Taking into consideration these different results we es- 
timate the diffusion of Vesta family members under the 
effect of Yarkovsky, YORP and collisions by a simple ran- 
dom walk process and the maximal diurnal Yarkovsky ef- 
fect between each reorientations. Fig. 9 shows a compari- 
son between the diffusion due to close encounters and the 
one due to the Yarkovsky effect as functions of the as- 
teroid diameter at different times (250 Myr, 1 Gyr and 3 
Gyr). For the close encounters, we plotted the diffusion 
obtained for the 11 asteroids and the extrapolation we de- 
duced for the whole main asteroid belt. For the Yarkovsky 
effect we plotted both the maximal diurnal effect (without 
any reorientation) and the result of the random walk pro- 
cess. We can see that reorienting collisions do not affect 
much the magnitude of the Yarkovsky effect for the range 
of diameters we are interested in, even at 3 Gyr. We re- 
call that the estimated age of the Vesta family is about 1 
Gyr (Carruba et al. 2007). For this duration, close encoun- 
ters are dominant for the diffusion of asteroids larger than 
40 ± 5 km (Fig. 9), but below this value the Yarkovsky 
effect is prevailing. For a diameter of 1 km the Yarkovsky 
effect is about 25 times greater than the effect of close 
encounters. 

Regarding the time evolution of these effects, we as- 
sume in this work that close encounters generate a ran- 
dom walk process that evolves as \/i whereas the maximal 
Yarkovsky effect is linear with time. When we take into 
account reorientations for the Yarkovsky effect we also 
have a random walk process so the evolution law should 
be proportional to \/t. However this is an asymptotic law 
which is valid when a great number of random walk steps 
is reached. In the case of reorienting collisions the char- 
acteristic time is about 300 Myr so there are only a few 
steps per Gyr and the asymptotic law is invalid on this 
time scale. It means that for the time scale we are inter- 
ested in, the Yarkovsky effect evolves faster than \/i so 
faster than close encounters. This is exactly what we ob- 
serve on Fig. 9 : for 250 Myr the close encounters take a 
more important part in the diffusion process (equivalence 
of the two effects for D=19 km) whereas for 3 Gyr the 
Yarkovsky effect is even more prevailing (equivalence for 
D=54 km). 
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Fig. 9. Comparison of the semi-major axis diffusion due to 
the close encounters and to the Yarkovsky effect after 250 
Myr (A), 1 Gyr (B) and 3 Gyr (C). EN is the diffusion due 
to close encounters with the f 1 asteroids considered in our 
simulations. TEN is the extrapolation we performed for 
the entire main belt (xl.3). YMX stands for the maximal 
diurnal Yarkovsky effect. YRW stands for the Yarkovsky 
effect with reorientations (random walk). 



6. Conclusion 

In this work we evaluate the effect of close encounters be- 
tween Vesta family members and massive asteroids on the 
semi-major axis diffusion of the family. We calculate this 
effect by two methods. We first look at the collective ef- 
fect of close encounters and deduce an estimation of the 
semi-major axis drift due to the 11 massive asteroids taken 
into account here. Then we study individually the close en- 
counters and make statistics over these events. With this 
second method we are able to separate the contribution 
of each of the 11 asteroids to the diffusion. We show that 
both methods give comparable results and that our find- 



ings are compatible with previous work on the Vesta fam- 
ily (Carruba et al. 2007). We use the results we obtained 
for the 11 asteroids to extrapolate and constrain the diffu- 
sion that would result from close encounters with all the 
main belt asteroids. We show that this diffusion should 
not exceed 1.3 times the diffusion due to the 11 asteroids. 
Finally, we compare the diffusion due to close encounters 
with the Yarkovsky-driven drift of the semi-major axis of 
asteroids in the 1-100 km range (in term of asteroids di- 
ameters). We find that both effects are equivalent over 1 
Gyr for a diameter of 40 ± 5 km. For smaller asteroids the 
Yarkovsky effect dominates the semi-major axis diffusion 
and for larger asteroids close encounters become more im- 
portant. Thus we confirm that although asteroids close en- 
counters have a significant influence, the main mechanism 
of semi-major axis diffusion that transports main belt as- 
teroids (and especially V-type asteroids) to Earth-crossing 
orbits via strong resonances is the Yarkovsky effect. 
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Abstract. We estimate numerically the semi-major axis chaotic diffusion of Vesta family asteroids induced by 

close encounters with 11 massive main belt asteroids : (1) Ceres, (2) Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) 
Hygiea, (15) Eunomia, (19) Fortima, (324) Bamberga, (532) Herculina, (704) Interanmia. Most of the diffusion is 
due to Ceres and Vesta. By extrapolating our results, we are able to constrain the global effect of close encounters 
with all the main belt asteroids. A comparison of this drift estimate with the one expected for the Yarkovsky 
effect shows that for asteroids whose radius is larger than about 35 km, close encounters dominate the Yarkovsky 
effect. Overall, our findings confirm the standard scenario for the history of the Vesta family. 



1. Introduction 

It is now admitted that most of V-type near-Earth aster- 
oids (NEAs) and howardite, eucrite and diogenite (HED) 
meteorites are fragments of a collision between Vesta and 
another object which also produced the Vesta family (see 
?). V-type NEAs and HED meteorites are former members 
of the Vesta family that have been carried from the main 
asteroid belt to Earth-crossing orbits by the two main res- 
onances in the neighborhood of Vesta : the and 3 : 1 
resonances with Jupiter. The major issue in this scenario 
is the fact that the average life time of fragments in these 
two resonances is too short to explain the the mean age 
of HED meteorites and V-type NEAs (?). 

The commonly accepted explanation for this is that 
NEAs spent most of their life time in the main asteroid 
belt before entering in resonance and being carried to their 
current orbits. This supposes that a mechanism can in- 
duce a diffusion process that brings fragments to one (or 
both) of the two main resonances. The Yarkovsky effect is 
now admitted to be the main mechanism explaining the 
diffusion of Vesta family members (e.g. ?). 

However other processes also contribute to the diffu- 
sion. It is in particular the case for the chaos induced 
by overlaps of mean motion resonances or close encoun- 
ters between asteroids in the main belt. In the case of 
the resonances, ? showed that the chaotic zone around 
the 3 : 1 resonance is too narrow to explain the needed 
diffusion. Recently, ? showed that close encounters among 
massive asteroids induce strong chaos which appears to be 
the main limiting factor for planetary ephemeris on tens 
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of Myr. In the continuation of this work we concentrated 
our study on the diffusion induced by close encounters of 
Vesta family members with the massive asteroids. 

This effect has already been studied for different as- 
teroid families (e.g. the Flora and Lixiaohua families ??). 
The most exhaustive study (in terms of number of aster- 
oids taken into account) concerned the Gefion and Adeona 
families (?). The author considered all 682 asteroids of ra- 
dius larger than 50 km for encounters with members of 
the families. They concluded that the four largest aster- 
oids ((1) Ceres, (2) Pallas, (4) Vesta and (10) Hygiea) had 
a much higher influence than all the 678 remaining aster- 
oids so the latters can be considered as negligible. The 
Vesta family has also been the object of such an analysis 
(?) but only close encounters with (4) Vesta were taken 
into account. These different studies show that the effect 
of close encounters with massive asteroids depends a lot 
on the considered family and its position in phase space 
with respect to those of massive asteroids. 

The purpose of this work is to explore further the 
effect of close encounters in the case of the Vesta fam- 
ily by considering more asteroids. We take into account 
the effect of the 11 following asteroids : (1) Ceres, (2) 
Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) Hygiea, (15) 
Eunomia, (19) Fortuna, (324) Bamberga, (532) Herculina, 
(704) Interamnia. Moreover, we use the results from these 
interactions to extrapolate and obtain an estimate of the 
effect which would be raised by the whole main asteroid 
belt. Finally we compare the semi-major axis diffusion in- 
duced by the Yarkovsky effect and the diffusion due to 
close encounters (in the case of the Vesta family) . 
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2. Numerical simulations 

We ran two sets of numerical simulations. Both of them 
comprise the 8 planets of the Solar System, Pluto, the 
Moon and the 11 asteroids enumerated above. In both sim- 
ulations, we generate a synthetic Vesta family produced by 
an initial collision of another object with Vesta. We do not 
model a realistic collision process but instead we consider 
a set of test particles initially at the same position as Vesta 
but with different relative velocities. This is not a problem 
because we are mainly interested in this work in exploring 
the phase space in the region of the Vesta family. The rel- 
ative velocities of the fragments of collision are sampled 
both in norm and direction. For the norm, we take 8 dif- 
ferent values (every 50 ms~^) between the escape velocity 
at Vesta's surface (around 350 ms~^) and twice this value 
(700 ms~^). The inclination is sampled between —50° and 
50° every 10° (11 different values). Finally, the direction 
in the invariant plane is taken every 10° (36 values). This 
way we create 3168 test particles that represent the Vesta 
family. We plotted the positions of these particles in the 
(a, e) and (o, /) plans (where a is the semi-major axis, e 
the eccentricity and / the inclination) superimposed with 
the catalog published by ? of current Vesta family mem- 
bers (Fig.l). Note that despite we did not compute the 
proper elements whereas the catalog uses them the two 
sets of points superimpose quite well. Moreover we con- 
structed an initial set of particles just after the original 
collision whereas the catalog exposes the remnants of the 
family after approximately 1.2 Gyr of evolution (see ?). 

In the first simulation (the reference simulation, S), 
the 11 asteroids are considered as test particles as well 
and thus do not perturb the fragments of collision. In the 
second simulation {SE) the 11 asteroids are considered 
as the planets and thereby interact with members of the 
synthetic Vesta family. Both solutions are integrated over 
30 Myr. 

!!!!!!!!!!!!!!!!!!!!!! Description mtegrateur !!!!!!!!!!!!!! 
In order to be able to analyze the evolution of the test 
particles and in particular the impact of close encoun- 
ters with the 11 asteroids we set up different logs. We 
record for each particle its instantaneous orbital elements 
every 1 Kyr. We do not compute proper elements but in- 
stead we use the minimum and maximum values (see ?) 
taken by a, e and / on time spans of 10 Kyr. Actually, we 
use the average of the minimum and maximum values as 
proper elements. We also set up a log of close encounters 
for each particle. Each time that a particle passes within 
0.01 AU from an asteroid, the asteroid number, the mini- 
mum distance of approach and the time of the encounter 
are recorded in the logs. 

3. Global overview of the diffusion 

The goal of this study is to understand whether the cur- 
rent flux of V-type asteroids coming from the main belt 
and carried to near-Earth orbits through strong planetary 
resonances can be explained by the chaos induced by close 
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Fig. 1. Initial distribution of the synthetic Vesta family 
in the semi-major axis - eccentricity plan (top) and the 
semi-major axis - inclination plan (bottom). Real mem- 
bers proper elements (from ?) are superimposed for com- 
parison. 



encounters in the main asteroid belt (or at least that these 
encounters contribute significantly). The simplest way to 
check this assumption is to look in both S and SE simula- 
tions to the number of NEAs as a function of time and to 
see if the consideration of asteroidal interactions implies 
a greater number of NEAs. We use the usual criterion to 
decide whether an asteroid is a NEA or not : if the perias- 
tron of its orbit become smaller than 1.3 AU the asteroid is 
considered as a NEA (see for example ?). In Fig. 2 we plot 
the evolution of the number of NEAs in both simulations 
as a function of time. We can see on this graph that the 
close encounters have not a major effect on the population 
of NEAs on the duration of the simulation. However, we 
did not ran our simulations on the total age of the Vesta 
family (around 1.2 Gyr) and we can think that just after 
the collision that originated the Vesta family, the popu- 
lation of V-type NEAs is dominated by fragments of the 
collision that were directly injected in the two main reso- 
nances. The diffusion process is supposed to provide these 
resonances on a much larger time scale. 

Another way to look for the effect of close encounters 
is to compare the proper elements of the initial and final 
distribution of fragments in both simulations. A usual way 
to do that is to compute the standard deviation of different 
proper elements (and in particular the semi-major axis) on 
the whole population of fragments at a time t with respect 
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Fig. 2. Number of near-Earth asteroids (NEA) among the 
synthetic Vesta family as a function of time in S and SE. 
The criterion for considering a fragment as a NEA is based 
on the minimum value taken by the periastron. When this 
value is less than 1.3 AU the fragment is considered as a 
NEA. 

to the initial values of these elements (see for example ?) 



N -I 



(1) 



where x can be the semi-major axis a, the eccentricity e 
or the inclination / of Vesta family fragments. As we men- 
tioned it we use the averages of the minima and maxima 
as proper elements and in particular for the initial values 
we compute them on the first 10 Kyr of the simulation. 
We do not use the initial conditions given in Fig.l because 
they are only instantaneous values. 

We want to measure the diS'usion in semi-major axis 
only due to close encounters, thereby we do not want to 
take into account resonances and in particular strong ones. 
For this purpose we plotted the value of the standard de- 
viation between the beginning and the end (after 30 Myr 
of evolution) of both simulations (S, SE) as a function of 
the initial (proper) semi-major axis. We divided the inter- 
val of semi major axis in bands of 0.01 AU and for each 
band we computed the standard deviation (of the semi- 
major axis, the eccentricity and the inclination) of the set 
of asteroids that are initially in this band. The results of 
this calculation are given in Fig. 3. We can see that the 
strong resonance 3 : 1 around 2.5 AU induces an impor- 
tant diffusion in both S and SE and that for semi-major 
axis lower than 2.25 AU both simulations are strongly af- 
fected by resonances (see ?, for a list of resonances acting 
on the Vesta family). Note that in strong resonances most 
fragments are highly unstable and for a great part of them 
the calculation stops before the end of the simulation (col- 
lisions with a planet or the Sun or escapes from the Solar 
System) . Thereby calculations of the diffusion between the 
beginning and the end of the simulation in such zones is 
not representative of their instability. For instance the dif- 
fusion is zero between 2.5 and 2.51 AU only because none 
of the fragments that were initially in this zone finished 
the simulation. 
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Fig. 3. Semi-major axis dependency of the diffusion in 
semi-major axis (top), eccentricity (middle) and inclina- 
tion (bottom) in S and SE. For each band of 0.01 AU, 
we plot the standard deviation of the average of the mini- 
mum and the maximum values of a (respectively e and I) 
during the last 10 Kyr of the simulations with respect to 
the initial values (see Eq.l). 



Between 2.26 and 2.48 AU, S shows a very limited 
semi major-axis diffusion even if we can see small peaks 
(e.g. around 2.42 AU corresponding to the 1 : 2 resonance 
with Mars). SE undergoes a more important diffusion in 
this band. We chose to run our calculations on this band 
because it fulfills at the same time two important condi- 
tions : the diffusion due to resonances is limited and the 
band contains (at the beginning of the simulation) a suf- 
ficient amount of fragments (2617 of the 3168 fragments) 
which is important for the statistics. Note that for the ec- 
centricity and the inclination, the diffusion is hidden in 
the remaining oscillations due to the secular evolution of 
these elements induced by the planets and that we can- 
not see any difference between both simulations in Fig. 3. 
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Fig. 4. Evolution in time of the semi-major axis diffusion 
in S and SE. We plot the standard deviation of the av- 
erage of the minimum and the maximum values of a for 
each fragment over 10 Kyr steps with respect to the ini- 
tial values (see Eq. 1) for asteroids whose semi-major axes 
are initially between 2.26 and 2.48 AU (in order to avoid 
strong resonances). 



However, this is not really a problem since in this work we 
are mainly interested in the semi-major axis diffusion. 

The evolution (in time) of the standard deviation of 
the semi-major axis computed on the 2.26 - 2.48 AU band 
for both simulations is plotted in Fig. 4. The standard de- 
viation is approximately constant for the reference sim- 
ulation which reinforces our choice of this band. On the 
contrary, the simulation which takes into account aster- 
oidal interactions clearly undergoes diffusion. Analyzing 
Fig. 4 gives us an approximation of the drift in semi-major 
axis due to asteroidal interactions during the simulation. 
After 30 Myr of evolution, we have a standard deviation 
around 9 x 10"* AU for SE and 2 x 10"* AU for S. In or- 
der to evaluate the diffusion due to asteroidal interactions, 
we have to remove the diffusion obtained in the reference 
simulation (which is probably due to resonances and out- 
put sampling of the secular evolution of the semi-major 
axis) from the total standard deviation. Note that the to- 
tal variance (the squared of the standard deviation) is the 
sum of the variances of the different contributions. Thus 
we obtain a standard deviation due to asteroidal interac- 
tions after 30 Myr of evolution of about 8.8 x 10"'' AU. 
We also observe several jumps in the semi-major axis dif- 
fusion, which are due to very close encounters of a single 
fragment with one of the 11 considered asteroids. Actually 
it has been shown (?) that the distribution of jump sizes 
is not a Gaussian but has thicker wings. Thus the curve 
is not regular on the time scale of this simulation because 
of the presence of rare very important events. We can find 
the same kind of effects on Fig. 3, the curve is not regular 
in the 2.26 - 2.48 AU band and shows peaks due to a few 
very close encounters. 
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Fig. 5. Semi-major axis evolution of a fragment which un- 
derwent various close encounters. We plot the maximum, 
minimum and instantaneous values of the semi-major axis 
of the fragment. All times of encounters closer than 10~^ 
AU are highlighted with black vertical lines. 
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Fig. 6. Semi-major axis evolution of a fragment which 
underwent a very close encounter (minimum distance of 
3.4 X IQ-^ AU) with Vesta (at t w 8.39 Myr). We plot 
the maximum, minimum and instantaneous values of the 
semi-major axis of the fragment. All times of encounters 
closer than 10~^ AU are highlighted with black vertical 
lines. 



4. DifTusion induced by close encounters among 
asteroids 

The global analysis that we made does not allow us to 
separate the contributions of the 11 asteroids in the diffu- 
sion process and to extrapolate to the whole main asteroid 
belt. In order to do that we have to examine logs of close 
encounters and estimate the diffusion resulting from each 
of these events and then make statistics. All this analysis 
is still restricted to the 2.26 - 2.48 AU band in order to 
avoid the strong resonances as previously. 

For each encounter recorded in the logs we compare 
the value of the semi-major axis before and after the en- 
counter. Fig. 5 gives an example of the evolution of the 
semi-major axis of a fragment that underwent several close 
encounters which resulted in jumps of different sizes. Fig. 6 
gives a more unusual example of the evolution of a frag- 
ment which underwent a very close encounter with Vesta 
resulting in a very important jump in semi-major axis. 
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Then for each of the 11 asteroids we calculate the stan- 
dard deviation of the distribution of jump sizes induced by 
close encounters with this asteroid. This gives a measure 
of the average diffusion resulting from a single encounter 
with the asteroid. As we already noticed, it is more con- 
venient to manipulate variances than standard deviations 
since the latter are not additive. If we assume that close 
encounters are independent events, every fragment under- 
goes a random walk and after TV encounters, the total 
variance is multiplied by a factor N : 

al[N]=Nal[l]. (2) 

This gives a measure of the average diffusion resulting 
from N encounters with the considered asteroid. We just 
have then to replace N by the mean number of close en- 
counters per fragment with the selected asteroid during 
the whole simulation and we obtain the diffusion (vari- 
ance) due to this asteroid on the duration of the simula- 
tion. The diffusion due to close encounters with all the 11 
asteroids on the duration of the simulation is the sum of 
the variances of the 11 asteroids. With the same reasoning 
than for Eq.2 it is possible to extrapolate the value of the 
variance on longer time scales : 

'^lit) = -^aliT^im). (3) 

sim 

Note that it as been shown (e.g. ?) that close encounters 
do not exactly result in a random walk (there are some 
correlations between successive encounters) and that the 
variance is not exactly linear with time (or number of 
encounters). It can, indeed, be written as a power law of 
the form : 

CTait) = Ct"". (4) 

For a random walk : B = 0.5. Numerical estimations 

of this coefficient B for different families (see ???) show 
that it depends on the family considered and can even 
evolve with time. However the values found in the litera- 
ture range between 0.5 and 0.65 and we did not find any 
estimations in the case of the Vesta family so we decided 
to stick to the random walk hypothesis. 

Note that in this reasoning we do not take into account 
the influence of the noise in the calculation of jumps sizes. 
Indeed, we observe on Fig. 5 and Fig. 6 that the minimum 
and maximum values of the semi-major axis are not con- 
stant outside the close encounters (jumps) but undergo os- 
cillations which are due to remaining secular evolution. It 
is possible to eliminate an important part of this noise by 
taking the minimum of minima and the maximum of max- 
ima over greater time scales before and after the encoun- 
ters. In order to choose the best time interval to compute 
these extrema and to evaluate the remaining noise after 
this treatment we use the reference simulation S and we 
calculate the standard deviation of the difference of semi- 
major £ixis before and after random times. Since there 
are no close encounters thus no jumps in this simulation 
this standard deviation measures only the remaining noise. 



With this method we compute the standard deviation of 
the noise for different time intervals for the calculation of 
extrema and we find that the minimal value (4.4 x 10~^ 
AU) is obtained for an interval length of 200 Kyr. This is 
the value we employ for all the following calculations. 

When we compute the variance of jump sizes in SE 
(for real close encounters) we are affected by the same 
noise. The value we obtain for the variance is the sum of 
the variance of the real diffusion and the variance of the 
noise. Thus the real diffusion resulting from close encoun- 
ters is given by : 

Another problem we experienced is that the logs of 
close encounters record all encounters within 0.01 AU but 
it appears that at this distance the effect is too small and 
is completely hidden in the noise. Thus statistics are con- 
taminated by false events (with no real jumps). In order 
to avoid this contamination we chose to take into account 
only encounters within 0.001 AU which are much more 
credible events. 

Finally, when a fragment undergoes several close en- 
counters spaced by less than 200 Kyr, we are not able to 
separate the effect of each encounter. We thus chose to ig- 
nore all such multiple encounters and we keep only single 
encounters in the calculation of the standard deviation of 
the size of the jumps. However, when we compute the total 
diffusion on the duration of the simulation we use the to- 
tal number of encounters that occur during the simulation 
as multiple encounters should not be ignored. 

Table 1 sums up the results we obtained for each aster- 
oid and for the 11 asteroids together. We give the numbers 
of encounter dtu-ing the simulation (including multiple en- 
counters), the standard deviation obtained from statistics 
on single encounters (cra[l])! the corrected values from 
noise (ffall]) and the standard deviations on the duration 
of the simulation {a-a{Tsim))- Finally, the percentages of 
contribution of the different asteroids are given in term of 
variance. We obtain a total diffusion due to close encoun- 
ters of 8.63 X 10"* AU during the 30 Myr of the simulation 
(see table 1). This number is to be compared with the drift 
rate obtained with the global approach (8.8 x 10~^ AU). 
The results given by both methods are in good agreement. 

By using Eq.3 we can extrapolate our results to longer 
times : 

aa{t) = 1.57 X IQ-^y/t (Myr) AU. (6) 

Of course it is possible to do the same for the contribution 
of each asteroid. This allow us to compare our results with 
the previous study of ?. In this article, the authors find a 
drift rate of 2.0l2:o ^ 10"^ AU/(100 Myr) by considering 
only close encounters with Vesta. If we extrapolate to 100 
Myr and consider only Vesta we find a drift of 1.3 x 10~^ 
AU which is coherent with ? findings. 

If we look at the different contributions in table 1 we 
see that close encounters with Ceres and Vesta together 
represent about 99% of the total diffusion (in variance). 
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Table 1 . Comparison of the contributions of the encoun- 
ters with the 11 considered asteroids to the semi-major 
axis diffusion. All standard deviations are given in AU. 
The contributions are given as percentages of the total 
variance (which is additive unlike the standard deviation) . 
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Contributions from the 9 other asteroids seem negligible. 
This is coherent with findings of previous studies on other 
asteroid families (see in particular ?). On the contrary, 
since Ceres represents 36% of the diffusion this asteroid 
should not be neglected (like in ?). We believe that two 
main reason can explain these proportions of contribu- 
tion : the mass of the asteroids and their proximity with 
the Vesta family members in phase space. The mass must 
influence the importance of the effect of a single close en- 
counter whereas the proximity in phase space influences 
the number of such events. In order to check for these 
correlations we plotted the variance due to close encoun- 
ters as a function of the mass of the considered aster- 
oid (Fig. 7) and the number of events as a function of the 
distance in phase space of the considered asteroid with 
respect to Vesta (Fig. 8). We use the definition of the dis- 
tance in phase space introduced by ? : 




-F2(5e)2-h2((5smi)2. 



(7) 



This distance has the dimensions of a velocity since it is an 
estimation of the ejection velocity that would be needed 
to carry two objects initially at the same position to their 
current positions. In Eq.7, we used the mean values of the 
semi-major axes, the eccentricities and the inclinations of 
the 11 asteroids during the whole simulation (30 Myr) in 
order to compute the distances. We deduce from these 
graphs that we indeed have an effect of the mass and the 
proximity in phase space but we cannot obtain a simple 
and precise law in order to evaluate the diffusion that 
would result from close encounters with other asteroids 
than the ones considered in our simulation. 

However, we can constrain the global effect of the re- 
maining objects of the main asteroid belt. In our simu- 
lation we took into account the most massive objects of 
the main asteroid belt. Actually, if we add the masses of 




Fig. 7. Contribution to the semi major axis diffusion due 
to each considered asteroid as a function of the mass of 
this asteroid. 




10 15 

d (AU/yr) 

Fig. 8. Number of close encounters (< 10~^ AU) between 
Vesta family fragments and each of the 11 considered as- 
teroids as a function of the distance in phase space be- 
tween these asteroids and Vesta. We use the mean values 
of the semi-major axes, the eccentricities and the inclina- 
tions of the 11 asteroids during the 30 Myr of the simula- 
tion in order to compute the distances (see Eq.7). 



the 11 asteroids of our simulation we obtain a total mass 
of 8.36 X lO~^*'Af0 whereas the estimated total mass con- 
tained in the main asteroid belt is about 15 x 10~^'^Mq 
(obtained from INPOPlOa fits, ?). The most massive as- 
teroid in the main belt is Ceres with 4.76 x IO-^^Mq, 
thus the remaining mass (not considered in our simula- 
tion) represents about 1.5 times the mass of Ceres. The 
total mass contained in the main asteroid belt is not well 
constrained and a previous determination of this mass by ? 
gave 18 X 10~^''Mq which implies a remaining mass closer 
to twice the mass of Ceres. Thus, it seems reasonable to 
fix the upper limit of this remaining mass to about twice 
the mass of Ceres. 

In Fig. 7, the slope of the regression curve is about 1.2. 
We plotted a straight line of slope 1 for comparison. A 
slope of 1 means that two objects of masses 0.5 have the 
same effect than one of mass 1. If the slope is higher than 
1 it is more efficient to have a single object whereas if 
the slope is lower than 1 it is more efficient to have two 
objects. Even if this slope is not well constrained in Fig. 7 
it is reasonable to consider that it is greater than 1 so 
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it is more efficient to have one big object tlian to spHt 
it into smaller objects. Tliereby, we can assume tfiat the 
total contribution of the remaining objects of the main 
belt will be less than twice the contribution of Ceres. 

We still have the problem of the influence of the dis- 
tance in phase space but it seems realistic to consider that 
the case of Vesta is singular since it is the parent body of 
the Vesta family and that the remaining objects must have 
a probability of close encounter closer to Ceres 's one than 
Vesta's one (see Fig. 8). Thus, the total variance due to 
encounters with all the asteroids of the main belt should 
not exceed (and is probably a lot smaller than) 1.7 times 
the value of the variance obtained in our simulation. In 
terms of standard deviations it corresponds to a factor of 
about 1.3. 



5. Comparison with the Yarkovsky-YORP effect 

The Yarkovsky effect is a non-gravitational force which re- 
sults from interaction between asteroids and the solar radi- 
ation. There are actually two versions of the Yarkovsky ef- 
fect : a diurnal and a seasonal effect. The diurnal effect can 
result in an increase or a decrease of the semi-major axis 
depending on the obliquity (e) of the asteroid (oc cose), 
whereas the seasonal effect systematically decreases the 
semi-major axis (cx sin^ e). The diurnal effect is maximal 
when e = 0° or 180° and is zero when e = 90°. On the 
contrary, the seasonal effect is maximal when e = 90° and 
zero when e = 0°, 180°. Both effects are size dependent, 
for instance the diurnal effect is given by (see in particular 
?) : 



The characteristic time scale of such events is given by 
(see ??) : 



$=2.5X10-^ 
dt 



1 km 
D 



cose AU/Myr. 



(8) 



In addition to the Yarkovsky effect, we have to consider 
the YORP (Yarkovsky - O'Keefe - Radzievskii - Paddack) 
effect which also comes from the solar radiation but affects 
the rotational velocity and the obliquity of asteroids. This 
effect depends a lot on the shape and the surface composi- 
tion of the asteroid but it is possible to do some statistics 
to understand what is the most probable scenario (see ??). 
It has been shown (?) that for most of basaltic asteroids 
(more than 95%), the obliquity is asymptotically driven 
to 0° or 180° by the YORP effect. The time scale for the 
YORP effect is about 10-50 Myr (see ??). This means that 
on a 1 Gyr time scale (which is the order of magnitude for 
the age of the Vesta family), we can consider that the 
obliquity is either 0° or 180° during the whole evolution 
and that the seasonal effect is negligible and the diurnal 
effect is maximal (in one or the other direction). 

Nevertheless, there is still another phenomenon which 
must be taken into account : the collisions. We have to 
distinguish whether a collision is destructive or not. The 
collisional life time of a 1-10 km asteroid in the main belt 
is about 1 Gyr. But non-destructive collisions are more 
frequent and can reorient the spin axis of the asteroid. 



Tcoll = 15 



R 

1 m 



1/2 



Myr = 335 



D 
1 km 



1/2 



Myr. (9) 



For kilometer-sized asteroids, the time between two reori- 
enting collisions is about 300 Myr, it means that there are 
a few (3 or 4) reorientations on the time scale of the age of 
the Vesta family. Thus it is likely that these reorientations 
will reduce the diffusion due to the Yarkovsky effect but 
not by large factor. Moreover, we can still consider that 
the YORP effect acts almost instantaneously between two 
reorientations and that the seasonal Yarkovsky effect is 
negligible and the diurnal effect is always meiximal. 

Taking into consideration these different results we es- 
timate the diffusion of Vesta family members under the 
effect of Yarkovsky, YORP and collisions by a simple ran- 
dom walk process and the maximal diurnal Yarkovsky ef- 
fect between each reorientation. Fig.9 shows a comparison 
between the diffusion due to close encounters and the one 
due to the Yarkovsky effect as functions of the asteroid 
diameter at different times of the evolution (250 Mjt, 1 
Gyr and 3 Gyr). For the close encounters, we plotted the 
diffusion obtained for the 11 asteroids and the extrapo- 
lation we deduced for the whole main asteroid belt. For 
the Yarkovsky effect we plotted both the maximal diurnal 
effect (without any reorientation) and the result of the 
random walk process. We can see that reorienting colli- 
sions do not affect much the magnitude of the Yarkovsky 
effect for the range of diameters we are interested in, even 
at 3 Gyr. We recall that the estimated age of the Vesta 
family is about 1 Gyr (?). For this duration, close encoun- 
ters are dominant for the diffusion of asteroids larger than 
35 km but below this value the Yarkovsky effect is prevail- 
ing. For a diameter of 1 km the Yarkovsky effect is about 
25 times greater than the effect of close encounters. 

Regarding the time evolution of these effects, we as- 
sume in this work that close encounters generate a ran- 
dom walk process that evolves as \/i whereas the maximal 
Yarkovsky effect is linear with time. When we take into 
account reorientations for the Yarkovsky effect we also 
have a random walk process so the evolution law should 
be proportional to \/t. However this is an asymptotic law 
which is valid when a great number of random walk steps 
is reached. In the case of reorienting collisions the char- 
acteristic time is about 300 Myr so there are only a few 
steps per Gyr and the asymptotic law is invalid on this 
time scale. It means that for the time scale we are inter- 
ested in, the Yarkovsky effect evolves faster than \/i so 
faster than close encounters. This is exactly what we ob- 
serve on Fig.9 : for 250 Myr the close encounters take a 
more important part in the diffusion process (equivalence 
of the two effects for D=19 km) whereas for 3 Gyr the 
Yarkovsky effect is even more prevailing (equivalence for 
D=54 km). 
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and that our findings are compatible with previous work 
on the Vesta family (?). We use the results we obtained 
for the 11 asteroids to extrapolate and constrain the diffu- 
sion that would result from close encounters with all the 
main belt asteroids. We show that this diffusion should 
not exceed 1.3 times the diffusion due to the 11 asteroids. 
Finally, we compare the diffusion due to close encounters 
with the Yarkovsky-driven drift of the semi-major axis of 
asteroids in the 1-100 km range (in term of asteroids diam- 
eters). We find that both effects are equivalent over 1 Gyr 
for a diameter of about 35 km. For smaller asteroids the 
Yarkovsky effect dominates the semi-major axis diffusion 
and for larger asteroids close encounters become more im- 
portant. Thus we confirm that although asteroids close en- 
counters have a significant influence, the main mechanism 
of semi-major axis diffusion that transports main belt as- 
teroids (and especially V-type asteroids) to Earth-crossing 
orbits via strong resonances is the Yarkovsky effect. 
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Fig. 9. Comparison of the semi-major axis diffusion due 
to the close encounters and to the Yarkovsky effect after 
250 Myr (top), 1 Gyr (middle) and 3 Gyr (bottom). EN is 
the diffusion due to close encounters with the 11 asteroids 
considered in our simulations. TEN is the extrapolation 
we performed for the entire main belt ( x 1.3). YMX stands 
for the maximal diurnal Yarkovsky effect. YRW stands for 
the Yarkovsky effect with reorientations (random walk). 



6. Conclusion 

In this work we evaluate the effect of close encounters 
between Vesta family members and massive asteroids on 
the semi-major axis diffusion of the family. We calculate 
this effect by two methods. We first look at the collec- 
tive effect of these close encounters and deduce an esti- 
mation of the semi-major axis drift due to the 11 massive 
asteroids taken into account here. Then we study individ- 
ually the close encounters and make statistics over these 
events. With this second method we are able to separate 
the contribution of each of the 11 asteroids to the diffu- 
sion. We show that both methods give comparable results 




